Transcriptome analysis of the transition from primary to secondary growth of vertical stem in Eucalyptus grandis.

Eucalyptus was one of the most cultivated hardwood species worldwide, with rapid growth, good wood properties and a wide range of adaptability. Eucalyptus stem undergoes primary growth (longitudinal growth) followed by secondary growth (radial growth), which produces biomass that is an important source of energy worldwide. In order to better understand the genetic regulation of secondary growth in Eucalyptus grandis, Transcriptome analyses in stem segments along a developmental gradient from the third internode to the eleventh internode of E. grandis that spanned primary to secondary growth were carried out. 5,149 genes that were differentially expressed during stem development were identified. Combining the trend analysis by the Mfuzz method and the module-trait correlation analysis by the Weighted Gene Co-expression Network Analysis method, a total of 70 differentially expressed genes (DEGs) selected from 868 DEGs with high connectivity were found to be closely correlated with secondary growth. Results revealed that the differential expression of these DEGs suggests that they may involve in the primary growth or secondary growth. AP1, YAB2 TFs and EXP genes are highly expressed in the IN3, whereas NAC, MYB TFs are likely to be important for secondary growth. These results will expand our understanding of the complex molecular and cellular events of secondary growth and provide a foundation for future studies on wood formation in Eucalyptus.

High-quality genome resource of Lasiodiplodia pseudotheobromae associated with die-back on Eucalyptus trees.

OBJECTIVES: Lasiodiplodia pseudotheobromae is an important fungal pathogen associated with die-back, canker and shoot blight in many plant hosts with a wide geographic distribution. The aim of our study was to provide high-quality genome assemblies and sequence annotation resources of L. pseudotheobromae, to facilitate future studies on the systematics, population genetics and genomics of the fungal pathogen L. pseudotheobromae. DATA DESCRIPTION: High-quality genomes of five L. pseudotheobromae isolates were sequenced based on Oxford Nanopore technology (ONT) and Illumina HiSeq sequencing platform. The total size of each assembly ranged from 43 Mb to 43.86 Mb and over 11,000 protein-coding genes were predicted from each genome. The proteins of predicted genes were annotated using multiple public databases, among the annotated protein-coding genes, more than 4,300 genes were predicted as potential virulence genes by the Pathogen Host Interactions (PHI) database. Moreover, the genome comparative analysis among L. pseudotheobromae and other closely related species revealed that 7,408 gene clusters were shared among them and 152 gene clusters unique to L. pseudotheobromae. This genome and associated datasets provided here will serve as a useful resource for further analyses of this fungal pathogen species.

Reference genes for Eucalyptus spp. under Beauveria bassiana inoculation and subsequently infestation by the galling wasp Leptocybe invasa.

Relative gene expression analysis through RT-qPCR is an important molecular technique that helps understanding different molecular mechanisms, such as the plant defense response to insect pests. However, the use of RT-qPCR for gene expression analysis can be affected by factors that directly affect the reliability of the results. Among these factors, the appropriate choice of reference genes is crucial and can strongly impact RT-qPCR relative gene expression analyses, highlighting the importance in correctly choosing the most suitable genes for the success of the analysis. Thus, this study aimed to select and validate reference genes for relative gene expression studies through RT-qPCR in hybrids of Eucalyptus tereticornis × Eucalyptus camaldulensis (drought tolerant and susceptible to Leptocybe invasa) under conditions of inoculation by the Beauveria bassiana fungus and subsequent infestation by L. invasa. The expression level and stability of eleven candidate genes were evaluated. Stability was analyzed using the RefFinder tool, which integrates the geNorm, NormFinder, BestKeeper, and Delta-Ct algorithms. The selected reference genes were validated through the expression analysis of the transcriptional factor EcDREB2 (dehydration-responsive element-binding protein 2). For all treatments evaluated, EcPTB, EcPP2A-1, and EcEUC12 were the best reference genes. The triplets EcPTB/EcEUC12/EcUBP6, EcPP2A-1/EcEUC12/EcPTB, EcIDH/EcSAND/Ecα-TUB, EcPP2A-1/Ecα-TUB/EcPTB, and EcPP2A-1/EcUPL7/EcSAND were the best reference genes for the control plants, mother plants, plants inoculated with B. bassiana, plants infested with L. invasa, and plants inoculated with B. bassiana and subsequently infested with L. invasa, respectively. The best determined reference genes were used to normalize the RT-qPCR expression data for each experimental condition evaluated. The results emphasize the importance of this type of study to ensure the reliability of relative gene expression analyses. Furthermore, the findings of this study can be used as a basis for future research, comprising gene expression analysis of different eucalyptus metabolic pathways.

Cytokinin promotes anthocyanin biosynthesis via regulating sugar accumulation and MYB113 expression in Eucalyptus.

Anthocyanins are flavonoid-like substances that play important roles in plants' adaptation to various environmental stresses. In this research, we discovered that cytokinin (CK) alone could effectively induce the anthocyanin biosynthesis in Eucalyptus and many other perennial woody plant species, but not in tobacco and Arabidopsis, suggesting a diverse role of CK in regulating anthocyanin biosynthesis in different species. Transcriptomic and metabolomic strategies were used to further clarify the specific role of CK in regulating anthocyanin biosynthesis in Eucalyptus. The results showed that 801 and 2241 genes were differentially regulated at 6 and 24 h, respectively, after CK treatment. Pathway analysis showed that most of the differentially expressed genes were categorized into pathways related to cellular metabolism or transport of metabolites, including amino acids and sugars. The metabolomic results well supported the transcriptome data, which showed that most of the differentially regulated metabolites were related to the metabolism of sugar, amino acids and flavonoids. Moreover, CK treatment significantly induced the accumulation of sucrose in the CK-treated leaves, while sugar starvation mimicked by either defoliation or shading treatment of the basal leaves significantly reduced the sugar increase of the CK-treated leaves and thus inhibited CK-induced anthocyanin biosynthesis. The results of in vitro experiment also suggested that CK-induced anthocyanin in Eucalyptus was sugar-dependent. Furthermore, we identified an early CK-responsive transcription factor MYB113 in Eucalyptus, the expression of which was significantly upregulated by CK treatment in Eucalyptus, but was inhibited in Arabidopsis. Importantly, the overexpression of EgrMYB113 in the Eucalyptus hairy roots was associated with significant anthocyanin accumulation and upregulation of most of the anthocyanin biosynthetic genes. In conclusion, our study demonstrates a key role of CK in the regulation of anthocyanin biosynthesis in Eucalyptus, providing a molecular basis for further understanding the regulatory mechanism and diversity of hormone-regulated anthocyanin biosynthesis in different plant species.

Identification of WUSCHEL-related homeobox gene and truncated small peptides in transformation efficiency improvement in Eucalyptus.

BACKGROUND: The WUSCHEL-related Homeobox (WOX) genes, which encode plant-specific homeobox (HB) transcription factors, play crucial roles in regulating plant growth and development. However, the functions of WOX genes are little known in Eucalyptus, one of the fastest-growing tree resources with considerable widespread cultivation worldwide. RESULTS: A total of nine WOX genes named EgWOX1-EgWOX9 were retrieved and designated from Eucalyptus grandis. From the three divided clades marked as Modern/WUS, Intermediate and Ancient, the largest group Modern/WUS (6 EgWOXs) contains a specific domain with 8 amino acids: TLQLFPLR. The collinearity, cis-regulatory elements, protein-protein interaction network and gene expression analysis reveal that the WUS proteins in E. grandis involve in regulating meristems development and regeneration. Furthermore, by externally adding of truncated peptides isolated from WUS specific domain, the transformation efficiency in E. urophylla × E. grandis DH32-29 was significant enhanced. The transcriptomics data further reveals that the use of small peptides activates metabolism pathways such as starch and sucrose metabolism, phenylpropanoid biosynthesis and flavonoid biosynthesis. CONCLUSIONS: Peptides isolated from WUS protein can be utilized to enhance the transformation efficiency in Eucalyptus, thereby contributing to the high-efficiency breeding of Eucalyptus.

Genes expression profiles in vascular cambium of Eucalyptus urophylla × Eucalyptus grandis at different ages.

BACKGROUND: Wood is a secondary xylem generated by vascular cambium. Vascular cambium activities mainly include cambium proliferation and vascular tissue formation through secondary growth, thereby producing new secondary phloem inward and secondary xylem outward and leading to continuous tree thickening and wood formation. Wood formation is a complex biological process, which is strictly regulated by multiple genes. Therefore, molecular level research on the vascular cambium of different tree ages can lead to the identification of both key and related genes involved in wood formation and further explain the molecular regulation mechanism of wood formation. RESULTS: In the present study, RNA-Seq and Pac-Bio Iso-Seq were used for profiling gene expression changes in Eucalyptus urophylla × Eucalyptus grandis (E. urograndis) vascular cambium at four different ages. A total of 59,770 non-redundant transcripts and 1892 differentially expressed genes (DEGs) were identified. The expression trends of the DEGs related to cell division and differentiation, cell wall biosynthesis, phytohormone, and transcription factors were analyzed. The DEGs encoding expansin, kinesin, cycline, PAL, GRP9, KNOX, C2C2-dof, REV, etc., were highly expressed in E. urograndis at three years old, leading to positive effects on growth and development. Moreover, some gene family members, such as NAC, MYB, HD-ZIP III, RPK, and RAP, play different regulatory roles in wood formation because of their sophisticated transcriptional network and function redundantly. CONCLUSIONS: These candidate genes are a potential resource to further study wood formation, especially in fast-growing and adaptable eucalyptus. The results may also serve as a basis for further research to unravel the molecular mechanism underlying wood formation.

Safety Assessment of the CP4 EPSPS and NPTII Proteins in Eucalyptus.

Glyphosate herbicide treatment is essential to sustainable Eucalyptus plantation management in Brazil. Eucalyptus is highly sensitive to glyphosate, and Suzano/FuturaGene has genetically modified eucalyptus to tolerate glyphosate, with the aim of both protecting eucalyptus trees from glyphosate application damage and improving weed management. This study presents the biosafety results of the glyphosate-tolerant eucalyptus event 751K032, which expresses the selection marker neomycin phosphotransferase II (NPTII) enzyme and CP4-EPSPS, a glyphosate-tolerant variant of plant 5-enolpyruvyl-shikimate-3-phosphate synthase enzyme. The transgenic genetically modified (GM) event 751K032 behaved in the plantations like conventional non-transgenic eucalyptus clone, FGN-K, and had no effects on arthropods and soil microorganisms. The engineered NPTII and CP4 EPSPS proteins were heat-labile, readily digestible, and according to the bioinformatics analyses, unlikely to cause an allergenic or toxic reaction in humans or animals. This assessment of the biosafety of the glyphosate-tolerant eucalyptus event 751K032 concludes that it is safe to be used for wood production.

Functional responses to multiple sequential abiotic stress (waterlogging-drought) in three woody taxa with different root systems and stress tolerance.

There is generally a trade-off in the resistance to drought and to waterlogging. However, several species are sequentially subjected to both stressors in many environments. We evaluated the ecophysiological strategies to cope with multiple sequential stress of waterlogging and drought (W + D) of three taxa differing in stress resistance and root morphology: the phreatophic Eucalyptus camaldulensis (Ec) and two shallow-rooted willow clones: Salix matsudana x Salix alba (SmxSa) and Salix nigra (Sn4). Individuals of the three taxa were grown in pots and assigned to either of four treatments: Control (well-watered plants), well-watered followed by drought (C + D); waterlogged for 15 days followed by drought (W15d + D) and waterlogged for 30 days followed by drought (W30d + D). Biomass allocation, growth (diameter, height, length of leaves, and roots), specific leaf area, stomatal conductance, water potential, hydraulic conductivity of roots and branches, leaf C13 and root cortical aerenchyma formation were determined at different stages of the experiment. Ec growth was not affected by W + D, developing tolerance strategies at leaf and whole plant levels. Differential effects of W + D were observed in both Salix clones depending on the time of waterlogging. In Sn4 and SmxSa, the root biomass was affected in W15d + D treatment, but a root tolerance response (aerenchyma and adventitious root formation) was observed in W30d + D. In the three taxa, and contrary to expectations, the previous exposure to a waterlogging period did not increase the susceptibility of the plants to a subsequent drought event. On the contrary, we found tolerance, which depended on the time of waterlogging exposure.

Genome-wide identification and expression profile analysis of metal tolerance protein gene family in Eucalyptus grandis under metal stresses.

Metal tolerance proteins (MTPs) as Me2+/H+(K+) antiporters participate in the transport of divalent cations, leading to heavy metal stress resistance and mineral utilization in plants. In the present study, to obtain better knowledge of the biological functions of the MTPs family, 20 potential EgMTPs genes were identified in Eucalyptus grandis and classified into seven groups belonging to three cation diffusion facilitator groups (Mn-CDFs, Zn/Fe-CDFs, and Zn-CDFs) and seven groups. EgMTP-encoded amino acids ranged from 315 to 884, and most of them contained 4-6 recognized transmembrane domains and were clearly prognosticated to localize into the cell vacuole. Almost all EgMTP genes experienced gene duplication events, in which some might be uniformly distributed in the genome. The numbers of cation efflux and the zinc transporter dimerization domain were highest in EgMTP proteins. The promoter regions of EgMTP genes have different cis-regulatory elements, indicating that the transcription rate of EgMTP genes can be a controlled response to different stimuli in multiple pathways. Our findings provide accurate perception on the role of the predicted miRNAs and the presence of SSR marker in the Eucalyptus genome and clarify their functions in metal tolerance regulation and marker-assisted selection, respectively. Gene expression profiling based on previous RNA-seq data indicates a probable function for EgMTP genes during development and responses to biotic stress. Additionally, the upregulation of EgMTP6, EgMTP5, and EgMTP11.1 to excess Cd2+ and Cu2+ exposure might be responsible for metal translocation from roots to leaves.

Survivorship and food consumption of immatures and adults of Apis mellifera and Scaptotrigona bipunctata exposed to genetically modified eucalyptus pollen.

Eucalyptus comprises the largest planted area of cultivated production forest in Brazil. Genetic modification (GM) of eucalyptus can provide additional characteristics for increasing productivity and protecting wood yield, as well as potentially altering fiber for a diversity of industrial uses. However, prior to releasing a new GM plant, risk assessments studies with non-target organisms must be undertaken. Bees are prominent biological models since they play an important role in varied ecosystems, including for Eucalyptus pollination. The main goal of this study was to evaluate whether a novel event (Eucalyptus 751K032), which carries the cp4-epsps gene that encodes the protein CP4-EPSPS and nptII gene that encodes the protein NPTII, might adversely affect honey bees (Apis mellifera) and stingless bees (Scaptotrigona bipunctata). The experiments were performed in southern Brazil, as follows: (i) larvae and adults were separately investigated, (ii) three or four different pollen diets were offered to bees, depending on larval or adult status, and (iii) two biological attributes, i.e., survivorship of larvae and adults and food intake by adults were evaluated. The diets were prepared with pollen from GM Eucalyptus 751K032; pollen from conventional Eucalyptus clone FGN-K, multifloral pollen or pure larval food. The insecticide dimethoate was used to evaluate the sensitivity of bees to toxic substances. Datasets were analyzed with Chi-square test, survival curves and repeated measures ANOVA. Results indicated no evidence of adverse effects of Eucalyptus pollen 751K032 on either honey bees or stingless bees assessed here. Therefore, the main findings suggest that the novel event may be considered harmless to these organisms since neither survivorship nor food consumption by bees were affected by it.

Transcriptome-Based Construction of the Gibberellin Metabolism and Signaling Pathways in Eucalyptus grandis × E. urophylla, and Functional Characterization of GA20ox and GA2ox in Regulating Plant Development and Abiotic Stress Adaptations.

Gibberellins (GAs) are the key regulators controlling plant growth, wood production and the stress responses in perennial woody plants. The role of GA in regulating the above-mentioned processes in Eucalyptus remain largely unclear. There is still a lack of systematic identification and functional characterization of GA-related genes in Eucalyptus. In this study, a total of 59,948 expressed genes were identified from the major vegetative tissues of the E. grandis × E. urophylla using transcriptome sequencing. Then, the key gene families in each step of GA biosynthesis, degradation and signaling were investigated and compared with those of Arabidopsis, rice, and Populus. The expression profile generated using Real-time quantitative PCR showed that most of these genes exhibited diverse expression patterns in different vegetative organs and in response to abiotic stresses. Furthermore, we selectively overexpressed EguGA20ox1, EguGA20ox2 and EguGA2ox1 in both Arabidopsis and Eucalyptus via Agrobacterium tumefaciens or A. rhizogenes-mediated transformation. Though both Arabidopsis EguGA20ox1- and EguGA20ox2-overexpressing (OE) lines exhibited better vegetative growth performance, they were more sensitive to abiotic stress, unlike EguGA2ox1-OE plants, which exhibited enhanced stress resistance. Moreover, overexpression of EguGA20ox in Eucalyptus roots caused significantly accelerated hairy root initiation and elongation and improved root xylem differentiation. Our study provided a comprehensive and systematic study of the genes of the GA metabolism and signaling and identified the role of GA20ox and GA2ox in regulating plant growth, stress tolerance, and xylem development in Eucalyptus; this could benefit molecular breeding for obtaining high-yield and stress-resistant Eucalyptus cultivars.

E. urophylla × E. grandis high-quality genome and comparative genomics provide insights on evolution and diversification of eucalyptus.

BACKGROUND: Eucalyptus urophylla × Eucalyptus grandis, an economically important forest tree, provides important raw material for energy and reduces damage to native forests. However, the absence of a high-quality E. urophylla × E. grandis reference genome has significantly hindered its evolution and genetic analysis. RESULTS: We successfully presented a high-quality reference genome of E. urophylla × E. grandis (545.75 Mb; scaffold N50, 51.62 Mb) using a combination of the Illumina, PacBio HiFi, and Hi-C sequencing platforms. A total of 34,502 genes and 58.56% of the repetitive sequences in this genome were annotated. Using genome evolution analyses, we identified a recent whole-genome duplication (WGD) event in E. urophylla × E. grandis. We further found that gene families associated with starch and sucrose metabolism, flavonoid biosynthesis, and plant-pathogen interaction were significantly expanded in E. urophylla × E. grandis. Moreover, comparative genomic and evolutionary analyses showed large structural variations among the different chromosomes of the 34 Eucalyptus accessions, which were divided into six clades. CONCLUSIONS: Overall, our findings provide a valuable resource for expanding our understanding of the E. urophylla × E. grandis genome evolution, genetic improvement, and its comparative biology.

High-throughput miRNA sequencing and identification of a novel ICE1-targeting miRNA in response to low temperature stress in Eucalyptus camaldulensis.

MicroRNAs (miRNAs) play a crucial role in the growth, development, morphogenesis, signal transduction, and stress response in plants. The ICE (Inducer of CBF expression)-CBF (C-repeat binding factor)-COR (Cold-regulated gene) regulatory cascade is an important signalling pathway in plant response to low temperature stress, and it remains unknown whether this pathway is regulated by miRNAs. In this study, high-throughput sequencing was employed for predicting and identifying the miRNAs that were likely to target the ICE-CBF-COR pathway in Eucalyptus camaldulensis. A novel ICE1-targeting miRNA, eca-novel-miR-259-5p (nov-miR259), was further analysed. A total of 392 conserved miRNAs and 97 novel miRNAs were predicted, including 80 differentially expressed miRNAs. Of these, 30 miRNAs were predicted to be associated with the ICE-CBF-COR pathway. The full-length of mature nov-miR259 was 22 bp and its precursor gene was 60 bp in length, with a typical hairpin structure. The RNA ligase-mediated 5' amplification of cDNA ends (5'-RLM-RACE) and Agrobacterium-mediated tobacco transient expression assays demonstrated that nov-miR259 could cleave EcaICE1 in vivo. Moreover, qRT-PCR and Pearson's correlation analysis further revealed that the expression levels of nov-miR259 were almost significantly negatively correlated with those of its target gene, EcaICE1, and the other genes in the ICE-CBF-COR pathway. We first identified the nov-miR259 as a novel ICE1-targeting miRNA, and the nov-miR259-ICE1 module may be involved in regulating the cold stress response in E. camaldulensis.

The In Planta Gene Expression of Austropuccinia psidii in Resistant and Susceptible Eucalyptus grandis.

Austropuccinia psidii, commonly known as myrtle rust, is an obligate, biotrophic rust pathogen that causes rust disease in a broad host range of Myrtaceae species. Eucalyptus grandis, a widely cultivated hardwood Myrtaceae species, is susceptible to A. psidii infection, with this pathogen threatening both their natural range and various forest plantations across the world. This study aimed to investigate the A. psidii transcriptomic responses in resistant and susceptible E. grandis at four time points. RNA-seq reads were mapped to the A. psidii reference genome to quantify expressed genes at 12 h postinoculation and 1, 2, and 5 days postinoculation (dpi). A total of eight hundred and ninety expressed genes were found, of which 43 were candidate effector protein genes. These included rust transferred protein 1 (RTP1), expressed in susceptible hosts at 5 dpi, and a hydrolase protein gene expressed in both resistant and susceptible hosts over time. Functional categorization of expressed genes revealed processes enriched in susceptible hosts, including malate metabolic and malate dehydrogenase activity, implicating oxalic acid in disease susceptibility. These results highlight putative virulence or pathogenicity mechanisms employed by A. psidii to cause disease, and they provide the first insight into the molecular responses of A. psidii in E. grandis over time.

Interspecies genome divergence is predominantly due to frequent small scale rearrangements in Eucalyptus.

Synteny, the ordering of sequences within homologous chromosomes, must be maintained within the genomes of sexually reproducing species for the sharing of alleles and production of viable, reproducing offspring. However, when the genomes of closely related species are compared, a loss of synteny is often observed. Unequal homologous recombination is the primary mechanism behind synteny loss, occurring more often in transposon rich regions, and resulting in the formation of chromosomal rearrangements. To examine patterns of synteny among three closely related, interbreeding, and wild Eucalyptus species, we assembled their genomes using long-read DNA sequencing and de novo assembly. We identify syntenic and rearranged regions between these genomes and estimate that ~48% of our genomes remain syntenic while ~36% is rearranged. We observed that rearrangements highly fragment microsynteny. Our results suggest that synteny between these species is primarily lost through small-scale rearrangements, not through sequence loss, gain, or sequence divergence. Further examination of identified rearrangements suggests that rearrangements may be altering the phenotypes of Eucalyptus species. Our study also underscores that the use of single reference genomes in genomic variation studies could lead to reference bias, especially given the scale at which we show potentially adaptive loci have highly diverged, deleted, duplicated and/or rearranged. This study provides an unbiased framework to look at potential speciation and adaptive loci among a rapidly radiating foundation species of woodland trees that are free from selective breeding seen in most crop species.

A Robust Disease Scoring Method to Screen Eucalyptus for Resistance Against the Aggressive Leaf Pathogen Teratosphaeria destructans.

Shoot and leaf blight caused by Teratosphaeria destructans is one of the most devastating foliar diseases on Eucalyptus. Therefore, breeding for resistance to this disease is considered urgent. Differences in susceptibility to T. destructans have been observed in the field but a robust inoculation protocol has, until recently, been unavailable and a disease scoring method for precise phenotyping has not been established. A first objective of this study was to determine the optimal conidial concentration for T. destructans inoculations on a susceptible Eucalyptus host. This concentration was then used to determine differences in susceptibility of six genotypes of Eucalyptus grandis × E. urophylla to the pathogen by assessing the percentage of infected stomata using electron microscopy and the percentage of leaf area covered by lesions (PLACL) using image processing. In addition, we developed a disease susceptibility index (SI) of six categories ranging from highly resistant (SI = 0) to highly susceptible (SI = 1.5 to 2). The more resistant genotypes were moderately resistant, with an SI value of 0.49 to 0.54 and a PLACL of 6.5 to 9%. In contrast, the more susceptible genotype scored an SI of 1.52 and PLACL of 48%. Host susceptibility was also assessed relative to the sporulation of the pathogen. This showed that the percentage of sporulation was not significantly correlated with host resistance. Overall, the results provide the basis for rigorous screening and selection of resistant genotypes to the disease caused by T. destructans using artificial inoculation.

Recruitment of distinct UDP-glycosyltransferase families demonstrates dynamic evolution of chemical defense within Eucalyptus L'Hér.

The economic and ecologically important genus Eucalyptus is rich in structurally diverse specialized metabolites. While some specialized metabolite classes are highly prevalent across the genus, the cyanogenic glucoside prunasin is only produced by c. 3% of species. To investigate the evolutionary mechanisms behind prunasin biosynthesis in Eucalyptus, we compared de novo assembled transcriptomes, together with online resources between cyanogenic and acyanogenic species. Identified genes were characterized in vivo and in vitro. Pathway characterization of cyanogenic Eucalyptus camphora and Eucalyptus yarraensis showed for the first time that the final glucosylation step from mandelonitrile to prunasin is catalyzed by a novel UDP-glucosyltransferase UGT87. This step is typically catalyzed by a member of the UGT85 family, including in Eucalyptus cladocalyx. The upstream conversion of phenylalanine to mandelonitrile is catalyzed by three cytochrome P450 (CYP) enzymes from the CYP79, CYP706, and CYP71 families, as previously shown. Analysis of acyanogenic Eucalyptus species revealed the loss of different ortholog prunasin biosynthetic genes. The recruitment of UGTs from different families for prunasin biosynthesis in Eucalyptus demonstrates important pathway heterogeneities and unprecedented dynamic pathway evolution of chemical defense within a single genus. Overall, this study provides relevant insights into the tremendous adaptability of these long-lived trees.

Haplotype mining panel for genetic dissection and breeding in Eucalyptus.

To improve our understanding of genetic mechanisms underlying complex traits in plants, a comprehensive analysis of gene variants is required. Eucalyptus is an important forest plantation genus that is highly outbred. Trait dissection and molecular breeding in eucalypts currently relies on biallelic single-nucleotide polymorphism (SNP) markers. These markers fail to capture the large amount of haplotype diversity in these species, and thus multi-allelic markers are required. We aimed to develop a gene-based haplotype mining panel for Eucalyptus species. We generated 17 999 oligonucleotide probe sets for targeted sequencing of selected regions of 6293 genes implicated in growth and wood properties, pest and disease resistance, and abiotic stress responses. We identified and phased 195 834 SNPs using a read-based phasing approach to reveal SNP-based haplotypes. A total of 8915 target regions (at 4637 gene loci) passed tests for Mendelian inheritance. We evaluated the haplotype panel in four Eucalyptus species (E. grandis, E. urophylla, E. dunnii and E. nitens) to determine its ability to capture diversity across eucalypt species. This revealed an average of 3.13-4.52 haplotypes per target region in each species, and 33.36% of the identified haplotypes were shared by at least two species. This haplotype mining panel will enable the analysis of haplotype diversity within and between species, and provide multi-allelic markers that can be used for genome-wide association studies and gene-based breeding approaches.

Natural annual transcriptome dynamics of Eucalyptus reveal seasonal adaptation of tropical/sub-tropical trees.

Seasonal environment cues are primary factors that influence a plant's growth and adaptation. The molecular basis of seasonal phenology has been well studied in trees growing in boreal and temperate ecosystems. However, little is known about the molecular phenology of trees belonging to tropical/sub-tropical regions. Here, we characterize the annual transcriptome dynamics of Eucalyptus dunnii, one of the world's most widely planted tropical/sub-tropical hardwoods, in natural environments. Our transcriptome analysis combined with the geographical distribution, environmental cues, microscopic observations and heterologous transformation analyses provides a molecular timetable of seasonal regulatory events of E. dunnii and its planting prospects in China. We further investigated the molecular mechanisms of the flowering phenology of E. dunnii. Our results suggest that low temperature is one of environment triggers for its seasonal flowering. In addition, a comparative transcriptome and cell ultrastructure analysis between Eucalyptus and Populus reveals the molecular bases of different shoot apex growth habits of trees originating from tropical/sub-tropical and boreal/temperate regions. Our study will provide cues for further investigating the molecular mechanisms underlying the seasonal phenology of trees from tropical/sub-tropical regions.

Comparative Transcriptome Analysis of High- and Low-Growth Genotypes of Eucalyptus urophylla in Response to Long-Term Nitrogen Deficiency.

Nutrients play important roles in the growth and development of most plant species. However, in perennial trees, the function of nutrients in different genotypes is poorly understood. Three different nutrient levels (low, sufficient, and high nutrient levels) were applied to two contrasting Eucalyptus urophylla cultivars (a high-growth cultivar ZQUA44 and a low-growth cultivar ZQUB15), and growth and expression levels were analyzed. Although the growth traits of both genotypes under nutrient starvation treatment were much lower than under abundant nutrients, tree height, crown width, and biomass of different ZQUA44 tissues were much higher than those of ZQUB15 at all three nutrient levels. Differentially expressed genes (DEGs) clustered into six subclusters based on their expression patterns, and functional annotation showed that the DEGs involved in glutathione metabolism and flavonoid biosynthesis may be responsible for nutrient starvation across different genotypes, while the DEGs involved in carotenoid biosynthesis and starch and sucrose metabolism may have a range of functions in different genotypes. The DEGs encoding the MYB-related family may be responsible for nutrient deficiency in all genotypes, while B3 may have different functions in different genotypes. Our results demonstrate that different genotypes may form different pathways to coordinate plant survival when they face abiotic stresses.